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Abstract Results on the lowest-order hadronic vacuum polarization contribution to the muon magnetic 
anomaly are presented. They are based on the latest published experimental data used as input to the 
dispersion integral. Thus recent results on r — > VrivK^ decays from Belle and on e^e~ annihilation to tt^tt" 
from BABAR and KLOE are included. The new data, together with improved isospin-breaking corrections for r 
decays, result into a much better consistency among the different results. A discrepancy between the Standard 
Model prediction and the direct g — 2 measurement is found at the level of 3(t. 
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1 Introduction 



The Standard Model (SM) prediction of the 
anomalous magnetic moment of the muon, a^i, is lim- 
ited in precision by contributions from hadronic vac- 
uum polarisation (HVP) loops. These contributions 
can be conveniently separated into a dominant low- 
est order {a}^'^'^'^) and higher order {a^^^'^-^'^) parts. 
The lowest order term can be calculated with a com- 
bination of experimental cross section data involv- 
ing e+e" annihilation to hadrons, and perturbative 
QCD. These are used to evaluate an energy-squared 
dispersion integral, ranging from the t:^^ threshold to 
infinity. The integration kernel strongly emphasises 
the low-energy part of the spectrum, dominated by 
the TTTT final stateEl When using e+e" data a devia- 
tion of more than Scr was observed [lil^ between the 
SM prediction and the direct experimental value 

A former lack of precise e"'"e~-annihilation data in- 
spired the search for an alternative. It was found 0] 
in form of t ^ i/^ + hadrons spectral functions, trans- 
ferred from the charged to the neutral state using 
isospin symmetry. During the last decade, new mea- 
surements of the TTTT spectral function in e+e~ an- 
nihilation with percent accuracy became available, 
superseding or complementing older and less precise 
data. With the increasing precision, which today 
is on a level with the t data in that channel, sys- 



tematic discrepancies in shape and normalisation of 
the spectral functions were observed between the two 
systems 0, Q- It was found that, when computing 
the hadronic VP contribution to the muon magnetic 
anomaly using the r instead of the e^e~ data for the 
27r and 47r channels, the observed deviation with the 
experimental value would reduce to less than Icr 
Fig. [1] summarizes the comparison between theory 
and experiment by 2006-8 
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Fig. 1. Comparison of the predictions for the 
muon magjnetic anomaly with the BNL mea- 
surement [4] in 2006, from (top to bottom) 

Refs. QlJEai- 

In this review I present the situation as of Octo- 
ber 2009, taking advantage of very recent papers: (1) 
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an updated analysis tllll u sing r data, including high- 
statistics Belle results [12] and an improved treatment 
of isospin-breaking corrections (IB) [l3|; (2) a BABMi 
measurement 1^ of the tttt spectral function using 



the hard initial state radiation (ISR) method, bene- 
fiting from a large cancellation of systematic effects in 
the ratio 7r7r7(7) to jj.ji'yipi) employed for the measure- 
ment; and (3) a global analysis [15| of all published 
e+e~ data. 



the use of the "dressed" cross section would entail 
the risk of double-counting some of the higher-order 
contributions. 

The function K{s) ~ 1/s in Eq. ^ gives a strong 
weight to the low-energy part of the integral. About 
91% of the total contribution to a'"'"*''"'-' is accumu- 
lated at centre-of-mass energies ^/s below 1.8 GeV 
and 73% is covered by the tttt final state, which is 
dominated by the p(770) resonance. 



2 HVP and g-2 



3 Updated 27r Analysis Using r Data 



It is convenient to separate the Standard Model 
(SM) prediction for the anomalous magnetic moment 
of the muon into different contributions, 



-,QED_ 



with 



,had,HO 



(1) 



(2) 



and where af-^ = (11 658471.810±0.016) 10"!° is the 
pure electromagnetic contribution |l6j . q}^'^'^'-' is the 
lowest-order HVP contribution, aJ;=>-d,HO ^ (-9.79 ± 
O.OScxp ± O.OSrad) 10^^° is the corresponding higher- 
order part and a):""'' = (15.4±0.1±0.2) lO-^^ 
where the first error is the hadronic uncertainty and 
the second is due to the Higgs mass range, accounts 
for corrections due to exchange of the weakly inter- 
acting bosons up to two loops [l3|- For the light- by- 
light (LBL) scattering part, (jh^-d.LBL^ ^^iq value 
(10.5±2.6) 10"^° from the latest evaluation 

Owing to unitarity and to the analyticity of the 
vacuum-polarization function, the lowest order HVP 
contribution to can be computed through the dis- 
persion integral [20| 



00 



(3) 



where K{s) is a well-known QED kernel, a — a{s = 0), 
and i?*"' (s) denotes the ratio of the "bare" cross sec- 
tion for e+e~ annihilation into hadrons to the point- 
like muon-pair cross section. The bare cross section 
is defined as the measured cross section corrected for 
initial-state radiation, electron-vertex loop contribu- 
tions and vacuum-polarization effects in the photon 
propagator. However, photon radiation in the final 
state is included in the bare cross section defined 
here. The reason for using the bare {i.e., lowest or- 
der) cross section is that a full treatment of higher 
orders is anyhow needed at the level of a^, so that 



3.1 Spectral functions in r decays 

The spectral function of the vector current decay 
T —i' X'l/r is related to the e~e~ — > X" cross section 
of the corresponding isovector final state X", 



f^xo (s) = wi,x-(s) , 



(4) 



where s is the centre-of-mass energy-squared or equiv- 
alently the invariant mass-squared of the r final state 
X, a is the electromagnetic fine structure constant, 
and Wi,x- is the non-strange, isospin-one vector spec- 
tral function corrected for IB and given by 



Vi,x~{s) 



6. 



1 dN. 



ds 



Nx 

-2 



(5) 



2s \ ^ Rm{s) 



with 



i?IB(s)^ 



FSR(5) Pl{s) 
Gem{s) PUs) 



F.{s) 



(6) 



In Eq. ([5]), {l/Nx)dNx /ds is the normalised invari- 
ant mass spectrum of the hadronic final state, and 
Bx- denotes the branching fraction of t — )• 
We use for the r mass the value rrir = (1776.84 ± 
0.17) MeV [21], and for the CKM matrix element 
]Kd] = 0.97418 ±0.00019 0, which assumes CKM 
unitarity. For the electron branching fraction we use 
B, = (17.818 ±0.032)%, obtained 23] supposing lep- 
ton universality. Short-distance electroweak radia- 
tive effects lead to the correction Sew — 1.0235 ± 
0.0003 @,[2i-[2i]. Ah the s-dependent IB corrections 
are included in Rib, which have been worked out in 
the 27r channel: FSR(s) refers to the final state ra- 
diative corrections 27|, and Gem{s) denotes the long- 



distance radiative corrections of order a to the pho- 
ton inclusive r spectrum, computing the virtual and 
real photonic corrections using chiral resonance [ist 
or vector dominance 1291. 
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3.2 Improved IB corrections 

As the physics of IB is described elsewhere [l^l I 
only concentrate here on the major difference between 
the previous analysis 0, 0] and the new one 11 1 . In 
Eq. ([6]) the ratio of the electromagnetic Fq{s) and 
weak F_ (s) form factors depends on the charged and 
neutral p parameters, as well as p — a; interference. 
While a small mass difference (of the order of 1 MeV) 
makes only a small effect on the dispersion integral 
because of its bipolar nature, a difference in the width 
can lead to a significant correction. Such a difference 
is expected from radiative p — tttt^ decays which have 
been evaluated in a scalar-QED vector-dominance 
model [s^l • The result is markedly different from the 
estimate made previously using only the hard radia- 
tion part [3l| . 

One could question the validity of using point-like 
pions in the calculation of radiative decays. How- 
ever several experimental tests support this assump- 
tion in e"'"e~ — > 7r"''7r~7(7) for the same mass range: 
lowest-order FSR with KLOE ^1, additional FSR 



with BaBAR [14 1 . 

The new IB corrections are listed in Table [T] 

3.3 Consistency of r spectral functions 

All published r 27r spectral functions are nor- 
malised to the world-average branching ratio. The 



shape of their mass dependence can be compared by 
looking at the relative difference between each spec- 
tral function and the combined one, locally averaging 
the data from ALEPH fisj, CLEO fsij, OPAL Q, 



and Belle [12[ (Fig.[2|). 



Table 1. Contributions to a),^'^'^'^ [7r7r,r] from 
the isospin-breaking corrections. Corrections 
shown in two separate columns correspond 
to the Gounaris-Sakurai (GS) and Kiihn- 
Santamaria (KS) form factor parametrisa- 
tions, respectively. 



Source 


A aj'^''^^ [7r7r,r] (IQ-^") 
GS Model KS Model 


Sew 


-12.21±0.15 


Gem 


-1.92±0.90 


FSR 


+4.67±0.47 


p — u) interference 


+2.80±0.19 +2.80±0.15 


m^±-m^o (cr) 


-7.88 


■m^± --m^o (Tp) 


-f4.09 +4.02 




+o.2o«:?^ +o.iit«:i? 


Trnj, EM decays 


-5.91±0.59 -6.39±0.64 



Total 


-16.07±0.59 


-16.70±0.64 
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Fig. 2. Relative comparison between the r — >■ UtTvtv^ spectral functions from ALEPH, CLEO, OPAL, Belle 
(data points) and their combined result (shaded band). 
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Since the world-average branching ratio is domi- 
nated by the ALEPH resuh, it is interesting to test 
the consistency between the absolute spectra, i.e. 
when each spectrum is normalised to the branching 
ratio measured by the same experiment. Fig. [3] shows 
a very good agreement between the full dispersion in- 
tegrals with comparable uncertainties. Thus the r 
experiments yield consistent absolute results. The 
average value and its error are rather insensitive to 
the branching ratio choice, although it is not true at 
the level of individual experiments. In particular the 
Belle result reaches its best precision only when the 
world-average is used. 



-c ALEPH 



-uCLEO 



tOPAL 



-u Belle 



I 



1 spectral functions 
combined 
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Fig. 4. Relative comparison between ee and r 
spectral functions, expressed in terms of the 
difference between neutral and charged pion 
form factors. Isospin-breaking corrections are 
applied to r data with the corresponding un- 
certainties included in the error band. 

As the g-2 dispersion relation involves an integral 
over the hadronic spectral function, it is interesting to 
consider the result with another kernel. By integrat- 
ing the e^e~ data weighted by the r matrix element 
C-r, and correcting for IB, one obtains the branching 
ratio B^^tF which can be directly compared to the 
measurements. Indeed, 



Fig. 3. Comparison of the aji^'^''^°[7r7r] values 
for different r experiments using their own 
r — > i/tTttt" branching ratios (closed circles) or 
the world-average (open circles). 



3.4 Comparison to e+e Data 

Figure |4] shows the relative difference between the 
ee and the IB-corrected t spectral functions versus s. 
The relative normalisation is consistent within the re- 
spective errors and the shape is found in better agree- 
ment than before 0], despite a remaining deviation 
above the p-mass-squared. The discrepancy with the 
KLOE data, although reduced, persists. 



dCVC 
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1 



4 m 2 

2s 

ml 



dss^^C^ (7) 



(8) 



The results for B^^o from e+e experiments 
are compared to the direct measurements in Fig. [S] 
The average, (24.78 ± 0.17exp ± 0.22ib)%, differ from 
the average t branching ratio, (25.42 ± 0.10)%, by 
(0.64±0.10^±0.17,^±0.22ji3)% to be compared to an 
applied IB correction of (-1-0.69)%. The discrepancy 
of about 2a is significantly reduced from the previous 
analysis [l[ (4.5(t). It should be emphasized that the 
observed increased deviation above the p mass be- 
tween ee and r spectral functions, essentially driven 
by the KLOE data, plays a more significant role in 
the B^^o'^ integral, rather than in the integral with 
its much steeper kernel. So we expect a better con- 
sistency for g — 2. 
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Fig. 5. The measured branching fractions for 
r — ^ VtIth^ (references in compared to the 
predictions from the e^e~ — >■ ■k'^-k~ spectral 
functions, applying the IB corrections. For the 
ee results, only the data from the indicated ex- 
periments in the 0.63 — 0.958 GeV range are 
used, and the combined ee data elsewhere. 
Vertical bands indicate average values 



4 Updated 2tx Analysis Using e+e In- 
cluding BaBAR 

4.1 The data 

Recent precision data, where all required radiative 
corrections have been applied by the experiments, 
stem from the CMD-2 fsTl and SND [si] experiments 
at the VEPP-2M coUider. They achieve comparable 
statistical errors, and energy-dependent systematic 
uncertainties down to 0.8% and 1.3%, respectively. 

These measurements have been complemented by 
results from KLOE [sgl at DA$NE running at the 
resonance centre-of-mass energy. KLOE applied for 
the first time the ISR technique to precisely determine 
the TTTT cross section between 0.592 and 0.975 GeV. 
The high statistics of the analysed data sample yields 
a 0.2% relative statistical error on the tttt contribution 
to a^^'^'^'~' . KLOE normalises the tttt^ cross section 
taking the absolute ISR radiator function from Monte 
Carlo simulation (Ref. 40] and references therein). 
The systematic error assigned to this correction varies 
between 0.5% and 0.9% (closer to the (/) peak). The 
total assigned systematic error lies between 0.8% and 
1.2%. 



In a recent publication IJ] the BaBAR Collabo- 
ration reported measurements of the processes ee — >■ 
7r7r7,/i/i7 using the ISR method at 10.6 GeV centre- 
of-mass energy. The detection of the hard ISR pho- 



ton allows BaBAR to cover a large energy range from 
threshold up to 3 GeV for the two processes. The 
7r7r(7) cross section is obtained from the 7r7r7(7) to 
MM7(7) ratio, so that the ISR radiation function can- 
cels, as well as additional ISR radiative effects. Since 
additional FSR photons are also detected, there is 
no additional uncertainty from radiative corrections 
at NLO level. Experimental systematic uncertainties 
are kept to 0.5% in the p peak region (0.6-0.9 GeV), 
increasing to 1% outside. 

4.2 Combining cross section data 

The details of the combination procedure are 
given in Ref. [l5| . The requirements for averaging 
and integrating cross section data are: (j) properly 
propagate all the uncertainties in the data to the fi- 
nal integral error, {ii) minimise biases, i.e., repro- 
duce the true integral as closely as possible in av- 
erage and measure the remaining systematic error, 
and (iii) minimise the integral error after averaging 
while respecting the two previous requirements. The 
first item practically requires the use of pseudo- Monte 
Carlo (MC) simulation, which needs to be a faithful 
representation of the measurement ensemble and to 
contain the full data treatment chain (interpolation, 
averaging, integration). The second item requires 
a flexible data interpolation method and a realistic 
truth model used to test the accuracy of the integral 
computation with pseudo-MC experiments. Finally, 
the third item requires optimal data averaging taking 
into account all known correlations to minimise the 
spread in the integral measured from the pscudo-MC 
sample. 

The combination and integration of the ee — )■ tttt 
cross section data is performed using the newly de- 
veloped software package HVPTools E6|- It trans- 
forms the bare cross section data and associated sta- 
tistical and systematic covariance matrices into fine- 
grained energy bins, taking into account to our best 
knowledge the correlations within each experiment as 
well as between the experiments (such as uncertain- 
ties in radiative corrections) . The covariance matrices 
are obtained by assuming common systematic error 
sources to be fully correlated. To these matrices are 
added statistical covariances, present for example in 
binned measurements as provided by KLOE, BABAR 
or the T data, which are subject to bin-to-bin mi- 
gration that has been unfolded by the experiments, 
thus introducing correlations. The interpolation be- 
tween adjacent measurements of a given experiment 
uses second-order polynomials, which is an improve- 
ment with respect to the previously applied trape- 
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zoidal rule. In the case of binned data, the interpo- 
lation function within a bin is renormalised to keep 
the integral in that bin invariant after the interpola- 
tion. The final interpolation function per experiment 
within its applicable energy domain is discretised into 
small (1 MeV) bins for the purpose of averaging and 
numerical integration. 

The averaging of the interpolated measurements 
from different experiments contributing to a given en- 
ergy bin is the most delicate step in the analysis chain. 
Correlations between measurements and experiments 
must be taken into account. Moreover, the exper- 
iments have different measurement densities or bin 
widths within a given energy interval and one must 
avoid that missing information in case of a lower mea- 
surement density is substituted by extrapolated infor- 
mation from the polynomial interpolation. To derive 
proper weights given to each experiment, wider av- 
eraging regions are defined to ensure that all locally 
available experiments contribute to the averaging re- 
gion, and that in case of binned measurements at 
least one full bin is contained in it. The averaging 
regions are used to compute weights for each exper- 
iment, which are applied in the bin-wise average of 
the original finely binned interpolation functions. If 
the value exceeds the number of degrees of free- 
dom (ndof), the error in the averaged bin is rescaled 
by V^JV^^^dof to account for inconsistencies. Fig. [6] 
shows the distributions in ^/s of the error recaling fac- 
tor, the relative weights for each experiment, and the 
contribution to the dispersion integral, as well as its 
error. It is seen that BaBAR dominates the averaging 
up to the p peak and above 0.95 GeV, while KLOE 
has a larger weight in-between owing to the steep be- 
haviour of the radiator function when approaching 1 
GeV. The uncertainty in the integral is dominated by 
the measurements below 0.8 GeV. 

The consistent propagation of all errors into the 
evaluation of a}^^'^-^'-' is ensured by generating large 
samples of pseudo experiments, representing the full 
list of available measurements and taking into ac- 
count all known correlations. For each generated 
set of pseudo measurements, the identical interpo- 
lation and averaging treatment leading to the com- 
putation of Eq. ([3]) as for real data is performed, 
hence resulting in a probability density distribution 
for ajj'''*''"°[7r7r], the mean and RMS of which define 
the IfT allowed interval. 



The fidelity of the full analysis chain (polyno- 
mial interpolation, averaging, integration) has been 
tested with toy models, using as truth representation 
a Gounaris-Sakurai vector-meson resonance model 
faithfully describing the tth data. Negligible biases 
below 0.1 (10^^° units) are found, increasing to 0.5 
(1.2 without the high-density BaBAR data) when us- 
ing the trapezoidal rule for interpolation instead of 
second order polynomials. 




Fig. 6. Top; Error rescaling factor accounting 
for inconsistencies among experiments versus 
^/s. Middle: Relative averaging weights for 
each experiment versus ^/s. Bottom: Contri- 
bution to the dispersion integral for the com- 
bined ee data obtained by multiplying the tttt 
cross section by the kernel function K{s) (solid 
line). The dashed (red) curve belonging to 
the right axis shows the corresponding error 
contribution, where statistical and systematic 
errors have been added in quadrature. 
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Fig. 7. Relative cross section comparison between individual experiments (symbols) and the HVPTools aver- 
age (shaded band) computed from all measurements considered. Shown are BABAR (top left), KLOE (top 
right), CMD-2 (bottom left) and SND (bottom right). 



The relative differences between BABAR KLOE, 
CMD-2, SND, and their average are given in Fig. [T] 
Fair agreement is observed, though with a tendency 
to larger (smaller) cross sections above ~0.8GeV for 
BaBAR (KLOE). These inconsistencies (among oth- 
ers) lead to the error rescaling shown in Fig. [6l 

4.3 Results for o^^'^'^'^Wtt] 

(J. L J 

A compilation of results for aj^'"*''"°[7r7r] for the var- 
ious sets of experiments and energy regions is given 
in Table [21 The inclusion of the new BaBAR data 
significantly increases the central value of the inte- 
gral, without however providing a large error reduc- 
tion, because of the incompatibility between mainly 
BABMi and KLOE, causing an increase of the com- 
bined error. In the energy interval between 0.63 and 
0.958 GeV, the discrepancy between the aj;'''*'^°[7r7r] 
evaluations from KLOE and BABAR amounts to 2. Oct. 
Since BaBAR is the only experiment covering the 
entire energy region between 2rn.^ and 1.8 GeV, it 
can provide its own evaluation JJJ of ajj''''''"°[7r7r] 
514.1±2.2,tat±3.1.y.t " 



more recent, but yet unpublished data [4l|. recovering 
agreement with the published SND cross sections 42 1 



Table 2. Evaluated aji^'^''^°[7r7r] contributions 
from the ee data for different energy intervals 
and experiments. Where two errors are given, 
the first is statistical and the second system- 
atic. The last value in parentheses is the total 
error). Also given is the r-based result. 







yfs (GeV) 


Exp. 


all^'^'-LOM (10-1°) 


2m^± — .3 


ee fit 


0.55±0.01 


.30 -.63 


Comb, ee 


132.6±0.8±1.0 (1.3) 


.63 -.958 


CMD2 03 


361.8±2.4±2.1 (3.2) 




CMD2 06 


360.2±1.8±2.8 (3.3) 




SND 06 


360.7±1.4±4.7 (4.9) 




KLOE 08 


356.8±0.4±3.1 (3.1) 




BABAR 


365.2±1.9±1.9 (2.7) 




Comb, ee 


360.8±0.9±1.8 (2.0) 


.958-1.8 


Comb, ee 


14.4±0.1±0.1 (0.2) 


Total 


Comb, ee 


508.4±1.3±2.6 (2.9) 


Total 


Comb, r 


515.2±2.0±2.7 (3.4) 



5 Multihadronic Contributions 



Since the new data are unavailable, we discard the 
obsolete CMD-2 data from the 7r7r27r° average, find- 
ing a;i-d,LO[^^27rO]=17.6±0.4etat±1.7,y,t (compared 



We also reevaluate the e+e — > tt+tt 2tt^ contri- to 17.0 ± 0.4stat ± 1.6syst when including the obsolete 



bution to a'^'''i'^°. It is found that the CMD-2 data 
used previously have been superseded by modified or 



CMD-2 data). The corresponding cross section mea- 
surements and HVPTools average are shown in Fig. [51 



tWhen not specified, the values are given in units of 10 
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From the still preliminary BaBAR results 43| it is clear 
that the region above 1.4 GeV is still underestimated 
at the present state, as corroborated by r data 0] . 



50 — . 



ND 



MSN 
- » DM2 
: ' OLYA 
40 — • SND 

Average 

30 — 




^Ji [GeVJ 

Fig. 8. Cross section measurements for 
e^e~ — > ■k'^-k~2-k^ used in the calculation of 
^had.LOj^^jvr"]. The shaded band depicts 
the HVPTools interpolated average within 
1(7 errors. The individual measurements are 
referenced in [gl. 



Results and comparison to experi- 
ment 



Adding to the ee-based aj;'"^''"'-'[7r7r] and 
^had,LO|-.^.^27rO] results the remaining exclusive multi- 
hadron channels as well as perturbative QCD [ij, we 
find for the complete lowest-order hadronic term 



a'^^'^-^o lee 



[ee] = 695.5±4.0,„p±0.7qcd (4.1tot). (9) 



It is noticeable that the error from the tttt channel 
now equals the one from all other contributions to 

^had.LO 

Adding further the other contributions (given in 
Section 2), we obtain the Standard Model prediction 
(stiU in 10"^" units) 

flSM [ee] = 11659 183.4±4.1±2.6±0.2 (4.9tot), (10) 

where the errors have been split into lowest and 
higher-order hadronic, and other contributions, re- 
spectively. The a^*^ [ee] value deviates from the ex- 
perimental average a"^^ = 116592O8.9ib5.4zb3.30, 
by 25.5 ±8.0 (3.2cr). For comparison the difference 
obtained with the updated r analysis is 15.7 ±8.2 
(1.9a). 

A compilation of recent SM predictions for 
compared with the experimental result is given in 
Fig. [HI The BABAR results are not yet contained in 



previous evaluations. The result by HMNT Q con- 
tains older KLOE data [i^, which have been super- 
seded by more recent results [39| leading to a slightly 
larger value for a„'"^''"°. 



HMNT 07 (e*e ) 

-276 ±51 I •— 

JN 09 (e'e) 

-290 ±65 I • 

Davier et al. 09 (i) 

-148 ±52 

Davier et al. 09 (e*e") 

-303 ±51 I • 1 

Davier et al. 09 (e*e" w/ BABAR) 

-246 ±49 I < 



BNL-E821 (WA) 
0±63 



«exp 



100 
X10-" 



Fig. 9. Compilation of recent results for a^^, 
subtracted by the central value of the experi- 
mental average Q. The shaded vertical band 
indicates the experimental error. The SM pre- 
dictions are taken from: HMNT 07 0], JN 
09 [i^, Davier et al. 09 (r-based and ee 
before BABAR) , and the ee-based value in- 
cluding BABAR. 



Conclusions: discussion and perpec- 
tives 

The following concluding remarks can be made: 

• The first point to emphasize is the better con- 
sistency between the ee and r analyses, result- 
ing from the improved IB corrections and the 
BaBAR results. There is still a difference of 
(6.8±2.9ee±3.4^) (1.5cr) in the tttt channel, but 
it can be considered as reasonable. The discrep- 
ancy was 2.9(7 before, reduced to 2.4cr after the 
T update. The other major difference affect- 
ing the two estimates is in the 27r27r'' channel, 
(3.8±1.7ee±l-4T) (1.7(t). For this case we have 
seen that a better measurement with e+e~ will 
likely move the ee result closer to the r value. 

• The internal consistency of the ee — ?■ tttt data 
is only fair, as a discrepancy is observed be- 
tween the BaBAR and KLOE results, which is 
increasing with energy. The difference on the p 
peak, about 3% (beyond the respective system- 
atic errors of 0.5% and 1.1%), is the most dam- 



-tThe g — 2 measurement is obtained from the ratio of two frequencies and needs as input the ratio of the muon to the proton 
magnetic moments. The latter ratio is derived from muonium hyperfine splitting, and its value has been updated 44] after the 
E-821 publication. The new value produces a shift of -1-0.92 10"^*^ of the Oi,. value. 
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aging for the dispersion integral. Here CMD-2 
agrees well with BaBAR, while SND lies between 
BaBAR and KLOE. However, the accuracy of 
both CMD-2 and SND is not enough to resolve 
the issue. 

• While the BaBAR measurement is explicitly 
done at NLO (including one additional ISR or 
FSR photon), it is insensitive to the Monte 
Carlo NLO generation. The situation is differ- 
ent for KLOE which relies on Phokhara for the 
ISR radiation function. Using the ISR process 
ee — >■ /i// BaBAR has been able to verify that 
Phokhara provides the right answer to an ac- 
curacy of 1.1%, however in a kinematic region 
(very hard ISR photons, x = 1 — s/sq > 0.9, 
where Sq is the square of the ee CM energy) far 
from that of KLOE (0.09 < a; < 0.66). A mea- 
surement of the muon ISR process by KLOE has 
been considered since some time and would be 
of considerable help to validate their approach. 

• Although the accuracy of the BaBAR results is 
similar to that of the combined previous ee ex- 
periments, the gain in precision that could have 
been hoped for was not realised because of the 
remaining discrepancy, essentially with KLOE. 

• All approaches now yield a deviation from the 
direct measurement, however at different lev- 
els depending on what mr input data are used: 
2.4ct with EaBAR alone, 3.2a with all ee data, 
3.7a with ee not including BaBAR, 2.9a with ee 
not including KLOE, and 1.9a with r alone. 

• Considering these results one can say there is 
some evidence for a deviation at the 3a level. 
The significance is still not enough to establish a 
breakdown of the Standard Model in the muon 
g—2, i.e. a contribution from new physics. How- 
ever it is a quantitative and very valuable in- 
formation which will help to constrain the new 
physics, if it is found at the LHC. 

This is the present situation. It should and will 
evolve as new results and new initiatives are taking 
place: 

• The hadronic spectral functions will continue 
to be refined as new results are expected from 
BaBAR in the multihadronic channels, from 
KLOE (tttt) with improved approaches (prelim- 
inary results are already available for the large- 
angle ISR analysis ^32] ) , and from the upgraded 



CMD-3 and SND detectors at the higher-energy 
VEPP-2000 collider [47|. 



• The theory error is still dominated by the un- 
certainty on the HVP contribution (4.1) , but 
it is now close to that of the hadronic LBL part 
(2.6)11. Since the latter contribution is unlikely 
to be known more precisely in the short-term it 
will eventually be the theory show-stopper. 

• But the real limitation at the moment is the g—2 
measurement itself. The uncertainty reached 
by E-821 is 6.3, larger than the full theory er- 
ror (4.9). It is therefore mandatory to pursue 
these measurements in order to reach higher 
precision. The factor 20 in precision obtained 
at BNL over the pioneering measurements per- 
formed at CERN has permitted to reach the 
electroweak scale in this process. Another fac- 
tor of 4, as anticipated by the new proposal 
submitted to Fermilab, or with the JPARC 



project [50|, will definitely provide quantitative 
information as we move into the new physics 
territory. 
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